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The low-energy singlet dynamics of the Quantum Heisenberg Antiferromagnet on the Kagome lattice is 
described by a quantitative Quantum Dimer Model. Using advanced numerical tools, the latter is shown to 
exhibit Valence Bond Crystal order with a large 3 6- site unit cell and hidden degeneracy between even and odd 
parities. Evidences are given that this groundstate lies in the vicinity of a Z 2 dimer liquid region separated by a 
Quantum Critical Point. Implications regarding numerical analysis and experiments are discussed. 
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The Kagome lattice, a two-dimensional corner- sharing ar- 
ray of triangles shown on Fig.[TJa), is believed to be one of the 
most frustrated lattices leading to finite entropy in the ground- 
state (GS) of the classical Heisenberg model Q. Hence, the 
Quantum S=l/2 Heisenberg Antiferromagnet (QHAF) on the 
Kagome lattice is often considered as the paradigm of quan- 
tum frustrated magnetism (2) where, in contrast to conven- 
tional broken symmetry phases of spin systems (as e.g. mag- 
netic phases), exotic quantum liquids or crystals could be re- 
alized. Among the latter, the algebraic (gapless) spin liquid is 
one of the most intriguing candidate Q. 

The herbertsmithite [4 ] compound is one of the few ex- 
perimental realizations of the S=l/2 Kagome QHAF. The ab- 
sence of magnetic ordering [ 5 ] down to temperatures much 
smaller than the typical energy scale of the exchange coupling 
J suggests that, indeed, intrinsic properties of the Kagome 
QHAF can be observed, even-though Nuclear Magnetic Res- 
onance and Electron Spin Resonance reveal a small fraction of 
non-magnetic impurities 1 6 ] and small Dzyaloshinsky-Moriya 
anisotropy 1 7 ] . So far, the nature of the non-magnetic phase is 
unknown and confrontation to new theoretical ideas have be- 
come necessary. Alternatively, ultra-cold atoms loaded on an 
optical lattice with tunable interactions might enable to also 
explore the physics of extended Kagome QHAF (8). 

The QHAF on the Kagome lattice has been addressed the- 
oretically by Lanczos Exact Diagonalization (LED) of small 
clusters (HHO). Despite the fact that the accessible cluster 
sizes remain very small, these data are consistent with a fi- 
nite spin gap and an exponential number of singlets within the 
gap, in agreement with a recent Density Matrix Renormaliza- 
tion Group study ifTTIl . In addition, an analysis of the four- 
spin correlations pointed towards a short-range dimer liquid 
phase lfT2l . Alternatively, a large-N approach fT3l and var- 
ious mappings to low-energy effective Hamiltonians within 
the singlet subspace |[T4lfT6l have suggested the formation of 
translation- symmetry breaking Valence Bond Crystals (VBC). 
Recently, recent series expansions around the dimer limit ifTTl 
showed that a 36-site VBC unit cell is preferred (see Fig[TJa)). 
In this context, the interpretation of the LED low-energy sin- 
glet spectrum remains problematic fT8ll . 

In this Letter, we use a quantitative dimer-projected effec- 




FIG. 1 : (color online) (a) Sketch of the 36-site VBC on the Kagome 
lattice showing "perfect hexagons" and (yellow) "star" resonances. 
Its unit cell is delimited by dashed (red) lines. The extension of 
the 108-site cluster (which fits exactly three 36-site cells) is (lightly) 
shaded and delimited by long-dashed (blue) lines. Hardcore dimers 
(formed by two neighboring spin-|) are located on the bonds (only 
shown on the shaded hexagons), (b) First Brillouin zone and avail- 
able momenta (labeling consistent with Ref. QjD). 



tive model to describe the low-energy singlet subspace of the 
Kagome QHAF. This is based on a controlled loop expan- 
sion fl9]| along the lines initiated by Rohksar and Kivelson 
(RK) [ 20 ] in the context of High-Temperature Superconduc- 
tivity and by Zeng and Elser I2H (ZE) for the Kagome lat- 
tice. LED of the effective model enable to reach unprece- 
dented sizes of periodic clusters which can accommodate a 
finite number of candidate VBC unit cells. The major results 
are (i) the numerical evidence of VBC order with a large 36- 



2 



site unit cell and a hidden degeneracy, (ii) the vicinity of a Z 2 
dimer liquid region separated by a Quantum Critical Point. 
Such results transposed to the original QHAF could explain 
some of its former puzzling numerical findings. 

Model - The method consists in projecting the QHAF, 
T-L = JJ2a j) Si • Sj, i and j being nearest neighbor (NN) 
lattice sites, into the manifold formed by NN Valence Bond 
(VB) coverings, an approximation shown to be excellent for 
the Kagome lattice [21 - 23 1. A transformation is then per- 
formed that turns the non-orthogonal VB basis into an orthog- 
onal Quantum Dimer basis. Overlaps between NN VB states 
can be written l24l (up to a sign) as a N ~ 2ni where N is the 
system size, n\ the number of loops obtained by superimpos- 
ing the two configurations, and a = l/y/2. This enables a 
systematic expansion in powers of a involving at order a 2p 
(p > 2) loops of sizes up to C = 2p + 2 (l9l|25) leading to a 
generalized Quantum Dimer Model (QDM) which, restricting 
to loops encircling only single hexagons, reads 



H (A[ = -Ja/o/ - Js (/~d7 + /Z^7 + ^§f) (l) 




where a sum over all the hexagons of the lattice of Fig. [TJa) 
is implicit. Kinetic terms (1) promote cyclic permutations of 
the dimers around the loops and diagonal terms (2) count the 
numbers of "flippable" loops. Here we use the approximate 
values (14-th order) , V 6 = 0.2, V 8 = 0.032, V 10 = V 12 = 0, 
J 6 = -0.8, J 8 = 0.251, Jio = -0.063, J 12 = (in units 
of J). Note that the kinetic processes Jc are quite close to 
ZE l2T1l initial lowest-order estimates, Jq = — 3/4, J 8 = 1/4 
and Jio = —1/16, a good sign of convergence of the expan- 
sion. Very importantly, we also include here the potential Vc 
terms which appear in our expansion scheme only at order 
2(C — 2) but play a major role. Note that an exact resumma- 
tion of the weights of (1), (2) up to all orders can be carried 
out (leading to tiny deviations) and that the "star" amplitudes 
J 12 and V12 vanish at all orders fT9l . 

Symmetry analysis - Although the generalized QDM ^ 
bears a sign problem, it can be addressed by LED of peri- 
odic 3 x L x L clusters which possess all the infinite lattice 
symmetries and L 2 unit cells. We shall consider L = 4 and 
L = 6 corresponding to the 48-site and 108-site clusters (see 
Fig.[TJa)), a tremendous improvement in terms of system size 
compared to the original QHAF, and use all available lattice 
symmetries (see available momenta in Fig. [TJb)) to block- 
diagonalize the Hamiltonian in its irreducible representations 
(IR) |26|. We further make use of a topological symmetry 
which splits the Hilbert space into 4 topological sectors (TS) 
{ni, 722, ^3}, where n a = (1) for an even (odd) number of 



dimers cut by each crystal axis a enclosing the torus l27l . 



VBC TAB CPE 

Af= 12 (+,+,±) (x,+,±) 

AT = 36 (+,+,±) (x,+,±) (+,x,±) (x,x,±) 

Af = 48 (+,+,±) (x,+,±) (x,x,±) (x,x,±) 

^2 (+,+,+) 



TABLE I: Quantum numbers of the GS multiplet for VBC with Jsf- 
site cells (JV = 3xraxmorA/' = 3x nV3 x n\^3) and degen- 
eracy J\f/3. The letters refer to the momenta of Fig. [T] Invariance 
under 27r/3 rotations (r3 = +) and parity under inversion {r^ = ±) 
and/or reflection about the momentum direction (a) are denoted as 
(7*3, 7*2, cr) consistently with Ref. fT8l . a = ± corresponds to even 
or odd GS and "x" means "symmetry not relevant". The Z2 dimer 
liquid (last line) GS in TS* is degenerate with the other TS GS l27l . 

An analysis of the low-energy spectrum and a careful in- 
spection of its quantum numbers provide invaluable informa- 
tions on the nature of the GS. For a VBC breaking discrete lat- 
tice symmetries, a multiplet structure is expected giving rise 
in the thermodynamic limit to a degenerate GS separated by a 
gap from the rest of the spectrum. We list in Table [I] the quan- 
tum numbers of the GS multiplet for the most popular VBC 
in the literature. 
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FIG. 2: (color online) (a) Complete low-energy excitation spectrum 
of the 48-site cluster as a function of the parameter A (see text). The 
lowest 25 levels are displayed in each TS (see (27)) using different 
symbols/colors. The 2-vison (blue lines) and 4-vison (dashed black 
line) gaps vanish around A = 1. (b) Close-up around A = 1 showing 
the excitation-energy of the lowest level of the most relevant IR (see 
Table[l| of the 108-site cluster: a QCP is identified from (i) a collapse 
of these excitations (b) and (ii) a very sharp peak in xa (c) (data for 
N = 108 divided by 200 to fit the scale). Note the transition is far 
more abrupt for N = 108 than for N = 48 (c). 

Quantum Critical Point - The QDM with only kinetic pro- 
cesses of equal amplitudes, Jq = J$ = Jio = Ji2(= 1/4) 
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provides an exactly solvable model 1291 Hrk with a gapped 
Z2 dimer liquid GS with short-range dimer correlations, sim- 
ilarly to the RK-point of the triangular QDM [30]. Interest- 
ingly, the first excitations (of energy J) correspond to pairs 
of (localized) topological vortices (visons). It is tempting to 
construct a simple interpolation H(\) = AH e ff + (1 — A)Hrk 
between this known limit and the effective model (1), (2). Its 
low-energy spectrum on the 48-site and 108-site clusters are 
shown as a function of A in Fig. [2ja) and Fig. |2jb), respec- 
tively. In the 48-site cluster, a high density of low-energy lev- 
els accumulate just above the GS at A ~ 1. This strongly 
suggests the vicinity of a Quantum Critical Point (QCP) char- 
acterized by vortex (vison) condensation as in the triangular 
QDM 13B . A closer look around A = 1 on the larger clus- 
ter reveals a sudden collapse of the low-energy excitations, 
clearly before A = 1. Remarkably, as shown in Fig.|2jc), this 
collapse coincides exactly with a very sharp peak of the sec- 
ond derivative of GS energy [28 ] xa = -d 2 (E% s /N)/d\ 2 , 
enabling to locate the QCP at A QC p ~ 0.9357. Note that the 
even-parity (cr = +) multi- vison excitations at momenta Ka 
and Kb merge with the GS precisely at Aqcp (see below). 
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FIG. 3: (color online) Complete low-energy excitation spectrum of 
the 48-site cluster as a function of J12. The lowest 25 levels are 
displayed in each TS (same symbols as Fig. |2ja)). The quantum 
numbers of the 4 lowest energy states for J12 < and J12 > are 
provided using different line types. 

VBC and hidden degeneracy - Next, to identify the ordered 
phase for A > Aqcp we add a finite star resonance amplitude 
J12 to i^ e ff. As seen in Figs. [3] and [4j this term has a crucial 
role. In fact, J 12 = is highly singular, with a (almost exact) 
degeneracy between odd (a = — ) and even (a = +) states. 
Physically, we believe it is related to the hidden Ising variables 
(introduced in Ref. 17) associated to the resonance parities of 
stars (see Fig.[TJa) for the 36-site VBC). Incidentally, it is re- 
markable that our effective model picks up such a feature via a 
vanishing effective J 12. For finite but still very small J12, this 



degeneracy is lifted (favoring a "ferromagnetic" Ising config- 
uration) allowing to characterize any candidate VBC (which 
should exist with both parity) from its lowest energy states. 

For the 48-site cluster, a close inspection of the associated 
quantum numbers and a comparison with Table [I] suggest a 
VBC with a 48-site (12-site) unit cell for J 12 < (J 12 > 0). 
For J 12 > slow fluctuations towards the pattern of a J\f = 48 
VBC are signaled by the presence of intermediate mid-gap 
states. Note that the lowest-energy excitations in the other TS 
(red lozenges) which barely depend on J 12 set roughly the 
(very small) VBC gap scale. Also, the high-density of levels 
within a small energy window of ~ 0.1 J above the GS is 
reminiscent of the singlet sector of the QHAF (T8J. 
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FIG. 4: (color online) (a) GS energies per site versus J 12 for 
K — Kr and a = =L The minimum defines the absolute GS energy 
E% s /N . Data for 48 sites (lines) and 108 sites (circles) are shown. 
An energy shift and scale transformation En — ^ \{En — N/2) al- 
low a direct comparison with twice the bond-energy of the QHAF. (b) 
Same as Fig.[2jb) but as a function of J 12 and compared to the N=48 
excitation at the T-point (full lines). Note the splittings between T, 
A and B levels of a few 10 -8 J are invisible at this scale. 



A similar analysis on the larger 108-site cluster l26l pro- 
vides definite evidence in favor of the 36-site (degenerate) 
VBC schematically depicted in Fig. []J a )> for A > Aqcp- 
First, Fig. |4ja) shows a kink at J 12 = of the GS energy 
per site, coinciding with the crossing of two GS of oppo- 
site parity and leading to almost the same slope discontinu- 
ity as for N = 48 (size effects are small). As for N = 48, 
this, in fact, corresponds to the level crossings of two groups 
of quasi- degenerate GS with opposite parities as shown in 
Fig. |4jb). We note an extremely fast decrease with N of the 
gaps between the A and T quasi-degenerate GS common to 
both clusters (a few 10" 3 J for = 48 to less than 10" 7 J 
for N = 108). Furthermore, the fact that a new state with 
momentum Kb (not allowed for N — 48) belongs to these 
two groups of quasi-degenerate GS definitely points toward a 
M = 36 unit cell (52 (which does not fit N = 48). This is 
further supported by the values of the average number Nc of 
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flippable length-^ loops compared in Table [II] to their values 
in VBC proposed in the literature. 



Groundstates 



N 6 /N H N 8 /N H N 10 /N H N 12 /N H 



N = 108; J12 = 0.01 0.154 0.274 0.491 0.081 

N = 108; J12 = -0.01 0.153 0.275 0.492 0.080 

3 X 2 X 2 |T4) 0.75 0.25 

3 X 2 V3 X 2^3 (EltBltni 0.167 0.25 0.5 0.083 



Z 2 dimer liquid |29 | 



1/32 15/32 15/32 1/32 



TABLE II: Average number Nc of flippable length-£ loops normal- 
ized to the total number Nh of hexagons. The last 3 lines correspond 
to the "frozen" limit of VBC and to the Z2 dimer-liquid. 

In summary, we introduced a generalized QDM to describe 
the low-energy physics of the QHAF on the Kagome lattice. 
In contrast to the latter, its groundstate properties can be ad- 
dressed by numerical simulations with unprecedented accu- 
racy for a frustrated quantum magnet. In particular, we pro- 
vide evidence of (i) a 3 6- site VBC order (with hidden degen- 
eracy), in agreement with recent series expansion ifTTlL and 
(ii) the vicinity of a QCP towards a topological Z2 dimer- 
liquid (cf. schematic phase diagram on Fig. 3 of Ref. IT9T ). 
Interestingly, a double Chern-Simons field-theory [33 ] also 
describes such a Quantum Critical Point. The above re- 
markable features of the generalized QDM transposed to the 
QHAF would resolve mysteries of the (small cluster) QHAF 
spectrum such as low-energy singlets carrying unexpected 
quantum numbers |18 ] and exceptional sensitivity to small 
perturbations (34). Experimentally, Kagome spin- 1/2 sys- 
tems generically contain small amounts of lattice and/or spin 
anisotropies (7) or even longer-range exchange interactions 
and the proximity to a QCP should render the experimental 
systems very sensitive to them. However, if under some con- 
ditions (pressure, chemical substitution,...) low-temperature 
spin-induced VBC order establishes, it could be revealed via 
small lattice modulations mediated by some magneto-elastic 
coupling. Lastly, we note that the above-mentioned perturba- 
tions as well as magnetic excitations [ 35 ] can be included in 
our scheme. 

We acknowledge support from the French National Re- 
search Agency (ANR) and IDRIS (Orsay, France). D.R 
thanks G. Misguich and A. Ralko for discussions. 



[1] J.T. Chalker, P.C.W. Holds worth, E.F. Shender, Phys. Rev. Lett. 
68, 855 (1992). 

[2] For review see e.g. "Two-dimensional quantum antiferromag- 

nets", G. Misguich & C. Lhuillier, in "Frustrated spin systems", 

Ed. H. T. Diep, World-Scientific (2005). 
[3] M. Hermele, Y. Ran, P. A. Lee, and X.-G. Wen, Phys. Rev. B 

77, 224413 (2008). 
[4] M. Shores, E. Nytko, B. Bartlett, and D. Nocera, J. Am. Chem. 

Soc. 127, 13462 (2005). 



[5] P. Mendels et al., Phys. Rev. Lett. 98, 077204 (2007). 

[6] A. Olariu et al., Phys. Rev. Lett. 100, 087202 (2008). 

[7] A. Zorko et al., Phys. Rev. Lett. 101, 026405 (2008). 

[8] J. Ruostekoski, Phys. Rev. Lett. 103, 080406 (2009). 

[9] P. Lecheminant et al., Phys. Rev. B 56, 2521 (1997). 
[10] P. Sindzingre et al., Phys. Rev. Lett. 84, 2953 (1999). 
[11] H. C. Jiang, Z. Y. Weng, and D. N. Sheng, Phys. Rev. Lett. 101, 

117203 (2009). 

[12] P.W Leung and Veit Elser, Phys. Rev. B 47, 5459 (1993). 
[13] J. B. Marston and C. Zeng, J. Appl. Phys. 69, 5962 (1991); 

SU(N) VBC have been introduced by N. Read and S. Sachdev, 

Phys. Rev. B 42, 4568 (1990). 
[14] A. V. Syromyatnikov and S. V. Maleyev, Phys. Rev. B 66, 

132408 (2002). 

[15] P. Nikolic and T. Senthil, Phys. Rev. B 68, 214415 (2003). 
[16] Ran Budnik and Assa Auerbach, Phys. Rev. Lett. 93, 187205 
(2004). 

[17] R.R.P. Singh and D.A. Huse, Phys. Rev. B 76, 180407 (2007). 
[18] G. Misguich and P. Sindzingre, J. Phys. Cond. Matt. 19, 145202 
(2007). 

[19] D. Schwandt, M. Mambrini andD. Poilblanc, arXiv: 1002.0774 
The weights of higher-order resonances (involving 2 or more 
hexagons) decay rapidly with C. 

[20] D.S. Rohksar and S.A. Kivelson, Phys. Rev. Lett. 61, 2376 
(1988). 

[21] Chen Zeng and Veit Elser, Phys. Rev. B 51, 8318 (1995). 
[22] M. Mambrini and F. Mila, Eur. Phys. J. B 17, 651 (2000). 
[23] For the frustrated QHAF on the square lattice see M. Mambrini, 

A. Lauchli, D. Poilblanc, and F. Mila, Phys. Rev. B 74, 144422 
(2006). 

[24] B. Sutherland, Phys. Rev. B 37, 3786 (1988); N. Read and 

B. Chakraborty, Phys. Rev. B 40, 7133 (1989). 

[25] For the same procedure on the frustrated QHAF on the square 
lattice see A. Ralko, M. Mambrini and D. Poilblanc, Phys. Rev. 
B 80, 184427 (2009). 

[26] E.g. the IR of momenta Kr (r3 — ri — +), Kb (^3 = 
+) and K A (r 2 = +) in TS* (27) of the 108-site cluster 
(a = ±) contain 79 548 096, 159 073 536 and 238 642176 
fully symmetrized states, respectively. Other momenta (with 
larger Hilbert spaces) have not been considered. However, apart 
from Kc, they should not be involved in the GS manifold (see 
Table|IJ. 

[27] Studied dimer liquids and VBC belong to the {0,0,0} 

[{1, 1, 1}] TS of the 48-site [108-site] cluster, named as TS*. 

"Other TS" refers to the 3 other (degenerate) TS. 
[28] As the fidelity, this quantity is a useful tool to identify quantum 

phase transitions. See S. Chen, L. Wang, Y. Hao and Y. Wang, 

Phys. Rev. A 77, 032111 (2008). 
[29] G. Misguich, D. Serban, and V. Pasquier, Phys. Rev. Lett. 89, 

137202 (2002). 

[30] R. Moessner and S. L. Sondhi, Phys. Rev. Lett. 86, 1881 (2001). 

[31] A. Ralko et al., Phys. Rev. B 76, 140404 (2007). 

[32] Reversely, the J\f = 48 cell does not fit in the 108-site cluster. 
The 3 x 4^/3 x 4\/3 cluster (144 sites) which can accommodate 
both J\f = 36 and M = 48 cells is out of reach. 

[33] C. Xu and S. Sachdev, Phys. Rev. B 79, 064405 (2009). 

[34] P. Sindzingre and C. Lhuillier, Eur. Phys. Lett. 88, 27009 
(2009). 

[35] A. Ralko, F. Becca and D. Poilblanc, Phys. Rev. Lett. 101, 

117204 (2008). 



